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Abstract 

In current lattice simulations of nucleon properties, the up and down quark 
masses are significantly larger than their physical values, while the strange 
quark can be included in simulations with its physical mass. When the up and 
down quark masses are much smaller than the strange-quark mass the chiral 
extrapolation of strange-quark matrix elements in the nucleon from the lattice 
up and down quark masses to their physical values can be performed with 
two-flavor chiral perturbation theory, thereby avoiding the slow convergence 
problem of the three-flavor chiral expansion. We explore the chiral expansion 
of several matrix elements of strange operators in the nucleon in two-flavor 
chiral perturbation theory and two-flavor partial-quenched chiral perturbation 
theory. 
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The strange quark, s, contribution to the properties of the nucleon provides a direct 
probe of the non- valence structure of the nucleon. If s-quarks were significantly more massive 
than Aqcd, as are the c- and 5-quarks, their contributions could be determined within the 
framework of perturbative QCD, and could be simply incorporated in the coefficients of 
operators involving up, u, and down, d, quarks and gluons, g. However, the s-quarks are 
light enough so that their contribution to nucleon properties are important and they should 
be considered as dynamical degrees of freedom in the nucleon, but heavy enough so that 
precise experimental measurements and the development of a rigorous theoretical description 
are challenging. 

Currently, several s-quark properties of the nucleon have been determined experimen- 
tally, including the contribution to the nucleon mass, m s (N\ss\N) Q, the s-quark helicity 
contribution to the proton spin, As 0, the s-quark contribution to the proton magnetic 
moment, fi s [|[, and the unpolarized s-quark parton distribution @]. Further, ongoing exper- 
iments are expected to map out the strange electric and magnetic form factors in the 
near future. For the determinations of m s (N\ss\N) and As, an expansion about the SU(3) 
flavor symmetry limit, in which the u-, d- and s-quarks are degenerate, is used but, unfortu- 
nately, in the absence of further experimental inputs, higher order SU(3) symmetry breaking 
effects |3| can only be estimated with hadronic models or with the leading non-analytic con- 
tributions from kaon loops computed in three-flavor chiral perturbation theory (xPT). To 
accommodate the sometimes slower than expected convergence of the SU(3) expansion in 
the baryon sector, conservative error estimates inevitably lead to large uncertainties. Lattice 
simulations of the s-quark properties of the nucleon, and the presently needed extrapolation 
to the physical quark masses must also deal with this issue. %PT [pj-Tl^] and its quenched- 



and partially-quenched extensions, quenched ;^PT (Q%PT) [0,14] and partially-quenched 
%PT (PQxPT) [|H)|-|H]], respectively, with both two and three flavors have been widely used 
to extrapolate physical quantities simulated on the lattice. One has confidence that for small 
enough quark masses both the two-flavor and three-flavor theories converge. However, in 
the three-flavor theory the relatively large value of the s-quark mass means that the issue of 
convergence must be explored on an observable-by-observable basis. If the SU(3) expansion 
is found not to converge rapidly it cannot provide a reliable extrapolation of the lattice 
calculations. However, the chiral extrapolations of strange matrix elements in the nucleon 
computed in lattice simulations do not suffer from the slow convergence of SU(3) %PT, and 
the two-flavor chiral expansion is sufficient for chiral extrapolations of these observables. 
The reason is very simple — one does not need to know the s-quark mass dependence for 
the purpose of chiral extrapolation. In lattice simulations only the u- and <i-quark masses 
are heavier than their physical values, while the s-quark can be simulated with its physical 
mass. Any variations in the s-quark mass can be implemented as an interpolation, rather 
than an extrapolation. 

In SU(2) ;\PT, the small expansion parameter of the theory is 

7TW A p 

e ~ A ' A ' A ' 1 } 

where in nature the physical pion mass is m n ~ 140 MeV, the A-nucleon mass difference 
is A ~ 300 MeV, and p is the characteristic momentum transfer in any given process. A 
is the "high-energy" scale set by the energy gap to strange hadronic excitations. A naive 
estimate of A ~ uik + rn^ — m N ~ 760 MeV, (accidentally close to the p meson mass) 
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roughly reflects the energy cost to create a ss pair. Thus e ~ 0.4 in the real world SU(2) 
expansion for low energy processes, consistent with the discussions found in Ref. [[12j| . In 
the chiral extrapolation of strange matrix elements, the e in eq. (ffl) also sets the radius of 
convergence of two-flavor PQxPT to be ~ 760 MeV. One expects the theory to converge 
rapidly for 'W< 300 MeV. 

For exploratory studies, several s-quark properties have been computed in quenched 
simulations fTEfl . As pointed out in Ref. the fact that large- N c quenched QCD (QQCD) 
(N c being the number of colors) and large- N c QCD are identical implies that differences in 
the low-energy constants of these two theories are 0(l/N c ). Thus after including quenched 
chiral logarithms, which are formally higher order in 1/N C but numerically enhanced by 
the smallness of light quark masses, one can hope that a QQCD calculation will have a 
0(1/N C ) ~ 30% error. However, for cases such as the proton strange magnetic moment, the 
higher order 1/N C corrections are likely to be important due to accidental cancellations at 
leading order (LO). In this work we present the chiral expansion of several nucleon strange 
matrix elements up to the next-to-next-to-leading order (NNLO) (i.e. 0(e 2 )) in both %PT 
and PQxPT [|i3|-|T7|j . PQxPT is presently the only tool that can rigorously remove (partial) 
quenching errors without performing fully unquenched simulations. We perform the PQ%PT 
calculations using two non-degenerate light-quark masses in the hope of providing an extra 
"handle" to aid in the extraction of strange matrix elements in the nucleon from lattice data. 



I. STRANGE MATRIX ELEMENTS IN QCD 

A. Spin Singlet Operators 

We start by considering the matrix element of the scalar current ss in x?T. This quark 
level operator is matched to the most general sum of hadronic operators in the heavy baryon 



chiral lagrangian JO]] which transform as spin and isospin singlets 

ss ^A S NN + B s T Q T a + • • • , (2) 

where N and T a are the nucleon and A-resonance fields, respectively, and where only the LO 
operators are shown. The prefactors can be rewritten as matrix elements of the quark oper- 
ator p0]| , and by taking the nucleon matrix elements at LO we obtain, A s = (N\ss\N)°, and 
B,s = — (A\ss\A)°, where the superscript denotes LO. At next-to-leading-order (NLO),there 
is a contribution proportional to A/A. The NNLO contribution comes from the one- loop 
diagrams shown in Fig. [I] which depend on both As and B$ operators, and tree-level coun- 
terterms. Thus up to NNLO (0(e 2 )), and working in the isospin limit, 

(N\ss\N) = (iVMAT - ((N\ss\N)° - (&\ss\A)°) J, 

where (Jajv = 1-8 is the AiW coupling, / ~ 132 MeV is the pion decay constant, m = m u = 
md is the quark mass in the isospin limit. The function J 7T =J(m n , A, /i) is 
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FIG. 1. One-loop graphs that give contributions of the form m q \ogm q to strange matrix el- 
ements in the nucleon. A single solid line denotes a nucleon, a double solid line denotes a 
/^.-resonance and a dashed line denotes a pion. Diagrams (a) and (b) are vertex corrections while 
diagrams (c) and (d) denote wavefunction renormalization. 



2 \ _ . r~r~n n, /A - VA 2 - m 2 + ie 



J(m, A, n) = (m 2 - 2A 2 ) log + 2AVA 2 - m 2 log ( 



A + VA 2 - m 2 



IE 



+ jj(m 2 -2A 2 ) , (4) 

where fi is the renormalization scale. The scale dependence of the counterterms Cs, Ds and 
Es exactly compensate the scale dependence introduced by the loop function, J n . 

An analogous result is found for the matrix elements of the s-quark spin-independent 
twist-2 operators 

where is the QCD covariant derivative and where (• ■ ■) indicates the symmetrization and 
removal of traces of the enclosed indices. We will first focus on forward matrix elements and 
leave the off-forward case to the next section. The operator O^ 1 '"* 1 ™ is matched to hadronic 
operators 

Op"'* 1 -> ai n) Nv (fil ■ ■ ■ v^N + (3i n) T a v ifM1 ■ ■ ■ v^T a 

+ 7 MT (M V 2 ■ ■ ■ V Vn~l T Vn) + . . . ? (6) 

where v M is the nucleon four- velocity, cui and the combination pi 1 — 7s /3 can be rewritten 
as the unpolarized nucleon and A twist-2 matrix elements, i.e. = (N\0^ 1 "' fln \N)° and 
p(n) _ 7 («)/ 3 = _(A|Of -^|A)°. The nucleon matrix element of Of is related to the 
n-th moment of the s-quark parton distribution 

(N\0^-" n \N) = 2(x n - 1 ) sN P^ 1 ---P>* n) , (n = 1,2,3,...) , (7) 
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where P M is the nucleoli momentum, and where 

(i""^ = f dx x n ~ l {s(x) + {-l) n s{x)) , (8) 



with s (s) is the quark (antiquark) spin-averaged distribution. The variable (Feynman) x is 
the momentum fraction of the proton carried by a quark in the infinite momentum frame, 
and for simplicity we have suppressed the renormalization scale dependence. To NNLO, 



A 

4tt 2 / 2 

+ c {n) — + D {n) — + E {n) — (9) 
s A s A 2 s A ' 

for n > 1. The n = 1 moment is not renormalized and remains zero, as it corresponds to 
the strangeness charge operator. It is no surprise that the chiral structure of (^ n_1 ) s ^ * s 

„n-l\ 



identical to the chiral structure of (x 11 ) u+dN that was found in Refs. |2l]j2!| and to that 
of (A^lssjA^) in eq.(|3|) as these are all isosinglet and spin-singlet operators. 



B. Spin Non-Singlet Operators 

The discussions of the previous section also apply to s-quark operators that are not spin 
singlets. We start with the proton strange magnetic moment induced by the strange vector 
current sj^s. This operator can be matched to tree level hadronic operators 

s^s - —e^v a d v N NSpN - 3/i° A T° S p T a ) + • • • , (10) 

771 j\r 

where Sp is the covariant spin operator, and where /x° N and /x° A are the LO strange mag- 
netic moments of the nucleon and A, respectively. At NNLO, (0(e 2 )), a straightforward 
calculation reveals 



o x 

H>s,N — ^ S ,N o 2 f 2 



A A 2 A v ' 

where g& ~ 1-26 is the ttN coupling constant, and L.,, = m 2 log (m 2 /^ 2 ). 
Consider now the matrix elements 

(JV l^^^siD* 42 • ■■iD^ ) s\ N) = 2 <x n_1 ) As N M N r {fll P^ ■ ■ ■ P^ , 
(N \sa a{ ^ l5 iD^ ■ ■ ■%D^s\ N) = 2 {x n ~ l ) &s N M N r [a P^ l] P» 2 ■ ■ ■ P^ , 

(n= 1,2,3,...) (12) 

where [• ■ •] means that the indices enclosed are antisymmetric and traceless. is the nucleon 
polarization vector with r 2 = —1. These matrix elements are related to moments of parton 
distributions 
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(^As 



dx X™" 1 (As (ar) + (-1) As (a;)) 



dx x 71 - 1 (Ss (x) + (-l) n 8s(x)) 



(13) 



where As (As) is the helicity distribution , and Ss (Ss) is the transversity distribution |2c 
of the s-quarks (s-antiquarks). (x n ~ 1 ) SsN and (x n ~ 1 ) AsN have the same chiral structure as 
[As, and up to NNLO 



(x n ~ l ) = (x n - l )° 

\ /As(8s),N \ /As{5s),N 



lTT 2 f 2 



A . A 2 



s{8s) 



(6s), I 



n A n ^A m 

+ <~>As(6s),N + ^As(5s),N + ^As(8s),N ~£ 



(14) 



However, it is important to note that not all spin non-singlet twist-2 matrix elements have 
this chiral structure. 

The s-quark contribution to the proton spin J s is related to the non-forward matrix 
element of the strange quark (traceless) energy momentum tensor 7^ [^|] which is just 
the operator shown in eq.([|) with n = 2, 



r {m,w) = s^HD^s . 



(15) 



This operator matches to |pq| 

^ 2^ 1 ^ a V<9, (J s j7V ~NSpN - 3J° A T^T,) + (x)° j7r d (w 7r a ^ 2) 7r a • ■ ■ , (16) 

in the hadronic theory, where ■/gjv(A) * s the s-quark contribution to the nucleon (A) 
spin, (x)° S7T is the LO momentum fraction of the pion carried by strange quarks and n a is 
the pion field (a = 1, 2, 3). It is straightforward to show that [^5| at NNLO, 



J. 



A 
A 



J «,N + C JsT + Dj S~\2 + Ej S 



Stt 2 / 2 



A 2 
A 2 

\0 



rn 
A 



(Js,n (^9 a L-k + 2g AN J n ) — (J[ 



s,A — o )yAN J TT 



• (17) 



II. STRANGE MATRIX ELEMENTS IN PQQCD 

In this section, we compute some of the quantities discussed in the previous section 



in PQxPT. The method for computation in PQxPT is well documented in Refs. 



especially in the case of two light, non-degenerate flavors [|17]]. In PQQCD, the fermion 
sector is generalized from two light quarks to four fermionic quarks (two "valence" 
quarks and two "sea" quarks) and two bosonic quarks ("ghosts"). The lagrangian has a 
SU(4\2)l Cg) SU(A\2)ji (g) U(l)y symmetry in the chiral limit and this symmetry is assumed 
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to be spontaneously broken down to S77(4|2)y <8> U(l)y so that a connection with QCD 
can be made. The valence quark masses and ghost quark masses are identical so that only 
the sea quarks contribute to the internal fermion loops (not coupled to electroweak gauge 
fields). If the ghost quark masses are set equal to infinity, the theory becomes the two-flavor 
quenched theory, but more importantly, if the sea quark masses are reduced to those of the 
valence quark masses, then QCD is recovered. In PQxPT, nucleons are embedded in a 70 
dimensional representation Bijk while the A-resonances are embedded in a 44 dimensional 
representation 7^, where the indices on both tensors run from k = 1 to 6. Therefore, 
at LO in PQ%PT the isoscalar operators that contribute to the matrix elements of interest 
have the form 

NYN^W k YB l3k , T a TT a ^T a ' t3k TT a , jk , (18) 

where r = 1 or S^. While it is true that the extension of the flavor structure of the strange 
operator from QCD to PQQCD is not unique, we have extended the flavor singlet under 
SU(2) to a singlet under SU(4|2) for simplicity. One could consider an extension to non- 
singlet operators, as discussed in Ref. |l6| , p6| , but in this case such an extension has no 
obvious advantages. 

From a practical point of view, the strategy we are suggesting is that simulations are 
performed with an unquenched s-quark, and two partially-quenched non-degenerate light 
quarks. This corresponds to PQQCD with SU(5|2) flavor symmetry. In order to perform 
the chiral extrapolation in the u and d quark masses one then matches the result of the 
SU(5|2) simulation onto SU(4|2) PQxPT. However, in order to make use of the hierarchy 
between the u, d and the s quark masses, both the valence and sea quark masses of the 
two light quarks must be much less than the mass of the s-quark. This approach should 
converge quite nicely for matrix elements in the nucleon for small enough meson masses, 
however, this is not expected to converge well for matrix elements in baryons containing one 
or more s-quarks. This is because the energy difference between an intermediate state with 
a kaon, in three-flavor xPT, e.g. £ — > NK — > S, and an intermediate state with a pion, e.g. 
£ — > S7r — ► S, is not large. 



A. Spin Singlet Operators 

The spin-singlet matrix elements we studied in the previous sections have the same chiral 
structures at NNLO. Since we have non-degenerate light quark masses, the s-quark matrix 
elements will be different for protons and neutrons at NNLO, but we will only consider the 
proton matrix elements here. We use the notation (0(2/+i,2S+i)) p to denote proton matrix 
elements of an operator with isospin I and spin S. Thus proton matrix elements of an 
isoscalar, spin-singlet operator, such as (^ n_1 ) sp and (p\ss\p) are denoted by . In 

PQxPT, up to NNLO, we find 



= (o M )° p - (<o (lll) >: - <e> (lll) >° A ) *b 



48vr 2 / 2 

+ C' SJ + D' s — + E' s ^ + Efo-j*. + E^^j^ , (19) 
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where rrij and mi are sea quark masses and the function Fq is 

F() ^Jud + Juu + Jju + Jlu + 2Jjd + 2Jid + 2'7r; tlj ^ u "I - 2T r j d ^ d ^ 

The functions J afe and 7; fe are J ab = J(m ab , A, fi) and 7; fe = H ab (J a , Jb, Jx), with 



n ab (A,B,C) = - 



{rrfj - m 2 Va ) (mg - m 2 ) (m| , - m 2 fc ) (mg - m 2 b ) 



4 

4 



A - 



(< - ) ( m i - m \ ) ( m t - m % ) - m x ) 
{m\ -m) 3 ){m\ - m 2 ) 

[m\ - m\ a ){m\ - m 2 J 



+ 



C 



B 



(20) 



(21) 



and the tree-level meson masses are mL, 



X (m x + m. 



and m^ 



m\). It is important to note that the low-energy constants Cg, D s and the Eg,i are related 



to those of xPT, as can be seen by considering the limit mj 
C s = C' s , D s = D' s , and E s = E' S I + E' s>2 + 2E' S>3 . 



m u and m; 



m d . i.e. 



B. Spin Non-Singlet Operators 

The observables fi s , (x n ~ l ) As and (x n ~ 1 ) Ss are matrix elements of operators of the form 
'C(i,3)) • Their chiral expansions, up to NNLO in PQ^PT, are 



^F + (/!^F 2 + ^ 2 F 3 + ^i^FoJ (0 (1 , 3) )° - ^g 2 AN F (0 (1)3) ) 



f22) 



127T 2 / 2 

where the loop functions are 



F\ = L u d + L uu + 2Lj U + 2L\ U + 3-R^^ 

F2 = 2L UU — L u d + Lj u + Li u + 3F^ U)J?ti + 3R nu;rid 

F3 = L uu — 5L u d + 2Li u + "SLid + 2Lj u + + 3R Vu:Vu + QR Vu:Vd + 3R nd:Vd , 

with L ab = m 2 b log (m 2 b //x 2 ), and R X} y = TC(L X , L y , Lx)- The counterterms Cy, Dy, and the 
Eyj are related to the appropriate counterterms in xPT in a straightforward way. Further, 
the s-quark contribution to the proton spin is 



J, 



s ,p 



12tt 2 / 2 



A 

A 



jo r i - n , A ' , F f m u , F > rn d m j + m t 

J s, P + C JsT ' D JsJ5 + h Js,i— + h Js,2— + h Js,r 



A 



A 



A 



g\F x + gi g A F 2 + ^g 2 F 3 + ^g 2 AN F 



5 , ( X )s,7T ^ 2 p 

'-^{ j s,a 2 — )9an^o 



(23) 



where the counterterms in eq. (|23| ) are related to the counterterms in %PT, given in eq. (IT 
by C Js = C' Js , D Js = D' Js , and E Js = f£ x + F^ + 2F (,) 



s 



III. CONCLUSIONS 



The chiral extrapolation of matrix elements computed in lattice QCD from the lattice 
quark masses to their physical values will play a central role in future lattice simulations. 
The slow convergence properties of three-flavor xPT will not impact the chiral extrapolation 
of strange matrix elements in the nucleon when simulations can be performed with pions of 
mass < 300 MeV. Two-flavor %PT is sufficient for the chiral extrapolation of these matrix 
elements. We have presented the chiral expansion of several strange matrix elements at 
NNLO in both xPT and PQxPT. 
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